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Abstract 

We study the statistics of the optimal path in both random and scale free networks, where 
weights w are taken from a general distribution P(w). We find that different types of disorder 
lead to the same universal behavior. Specifically, we find that a single parameter (S = AL~^^'^ for 
d-dimensional lattices, and S = AN^^^^ for random networks) determines the distributions of the 
optimal path length, including both strong and weak disorder regimes. Here z/ is the percolation 
connectivity exponent, and A depends on the percolation threshold and P{w). For P{w) uniform, 
Poisson or Gaussian the crossover from weak to strong does not occur, and only weak disorder 
exists. 
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The study of the optimal path in disordered networks has attracted much interest in 
recent years, due to its relation to many complex systems of interest, including polymers 
0], nanomagnet transport P], surface growth Q], spin glasses 0, tumoral growth immune 
response jol] , and complex networks [7] . The disorder is realized by assigning random weights 
w to the links of the network, where the weights are drawn from a distribution P{w), 
chosen to reflect the characteristics of the system under study. Examples include exponential 
disorder for tu„,re..ing effects g Power law d.sorder for driven diflusiorr in random media g 
lognormal distribution in conductance of quantum dots M and Gaussian distribution for 

n 

polymers [10|. 

Of particular current interest is the scaling behavior of the optimal path length i for 



different weight distributions 



The optimal path between two nodes is the path that 
minimizes the sum of the weights along the path. By definition, the weak disorder limit 
is obtained when almost all links contribute to the total weight of the optimal path. For 
this case ^ ~ L [in| for a ci-dimensional lattice of linear size L, and i ~ logA^ for random 



networks of N nodes |2I|- The strong disorder limit is obtained when a single link weight 
dominates the sum of weights along the path[ll|. For this case 



^L'^op', (1) 

with dopt = 1-22 Jd = 2) and dopt = 1-42 {d = 3) for lattices , and £ ~ AtVs 

random networks 7|. It is commonly agreed that strong disorder arises only when P(w) is 
broad enough. 

Lacking is a general criterion to determine which form of P{w) can lead to strong disorder, 
and a general condition when strong or weak disorder occurs. Here we derive analytically 
such a criterion, and support our derivation by extensive simulations. Using this criterion 
we show that certain power law distributions and lognormal distributions P(w) can lead to 

n 

strong disorder and to a weak-strong disorder crossover We also show that for P{w) 

uniform, Poisson or Gaussian, only weak disorder occurs regardless of the broadness of P{w). 
Importantly, we find that for all P{w) that possess a strong- weak disorder crossover, the 
distributions p{i) of the optimal path lengths display the same universal behavior. 

A useful way to study properties of the P{w) is to relate it to the uniform distribution 
tt{x) = 1, with X G [0, 1). From the relation P{w)dw = TT{x)dx = dx, w can be obtained 
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through the "disorder function" f{x) defined as w = f{x), which is equivalent to 

X = f-\w) = / Piw')dw'. (2) 
Jo 

In this form, every hnk has assigned to it a number x between and 1, and its corresponding 
weight is given by the relation w = f{x). Using /(x), we can write the total weight Wopt = 
^j=i'Wj of the optimal path as 

Wopt = fixi) + f{x2) + ... + fixi), (3) 

where the terms of the sum are arranged in decreasing order of values of wj = f{xj)- 

Next, we find the condition for strong disorder. We consider the first two terms of Eq. Q, 
since for strong disorder only the first term of is dominant ^ . The reason for choosing 
only two terms is that if the second term is very small compared to first then all others can 
be neglected and we are in the strong disorder, while if the second term cannot be neglected 
all others also cannot be neglected and we are in the weak disorder. Thus only the first two 
terms determine the disorder. By definition W2/W1 = f{x2)/f{xi) < 1, so we can expand 
/(X2) around Xi, 



Wi+W2 = f{xi) 

where 

S = 



1 , 

d{\nf) 



/(^i)[2-^], (4) 



dx 



{X1-X2). (5) 

X=Xl 

To see that S determines the disorder strength and the weak-strong disorder crossover, we 



note that (i) for the strong disorder limit, Wopt ~ fi^i), and therefore f{x2)/f{xi) 0. 
Since S" ~ 1 — /(x2)//(xi), it follows that S" — 1, and (ii) for the weak disorder limit all 
weights Wj contribute to Wopt, so /(xi) > f{x2), and 5* <^ 1. 

To identify Xi and X2 in (0), we show that the strong disorder case is related to percola- 
tion ^12^. The percolation problem corresponds to the construction of random networks by 
the removal of links with probability 1— p. The removal of links relates to the strong disorder 
limit, because all links with weights Wj larger than necessary to keep the network connected 
are not relevant and can be removed [l3|. Therefore, the weights that are typically removed 
are those satisfying w > f{pc), where pc is the critical percolation threshold, the value of p 
for which the system switches between being connected to disconnected. Hence the value of 
Xi should be pc (in the limit of infinite system size). 
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For percolation on a finite system, pc is not unique but follows a distribution characterized 
by a width which scales as L~^/^ , where v is the percolation connectivity length exponent 



15|. Thus for strong disorder, both Xi and X2 are distributed close to Pc with widths L . 



For typical values of X\ and X2, x\ and x\ respectively, we expect 



~ L-^l\ (6) 

Vc 

Combining Eqs. ^ and (jH)), we find 

rf(ln/) 

Using w = /(x), and f'{x) = 1/ P{w)\w=f{x), Eq- O can be expressed as 

S ^ ^^^^ ^ = (8) 
where i^c is the solution of the equation 

^ P{w)dw = Pc. (9) 



Equation (jSJ expresses 5* as a function of the general disorder distribution P{w) and L. 

Next we extend our theory for S to random networks. Percolation at criticality on 
Erdos-Renyi (ER) networks [l6| is equivalent to percolation on a lattice at the upper critical 
dimension 4 = 6 Q- For = 6, L ~ A^^^ and u = 1/2. Thus, using (jHI) we obtain, 

S ~ AL'^/'' = AN-^/'^" = AN-^/\ (10) 

Thus, for ER networks we simply replace L~^^'^ in (jSJ with N^^^'^. 

Following similar arguments for a scale-free network with degree distribution g(k) ~ /c^''^ 
and 3 < A < 4, we can replace L-^/" by A^-(^-3)/(a-i) gjnce 4 = 2(A-l)/(A-3) [17]. Thus, 

S ~ AL-^^"" = AN-^/^'^^''^ = AAr-(^-=^)/(^-i). (11) 

For A > 4, Eq. (HH) reduces to (HU)) Q- 

Our main results, Eqs. (jH)), (fTUj) and (fTTj). suggest that since S controls the disorder, for 
the same S, we expect to obtain the same optimal path behavior independent of P{w). 

Next we calculate A for several specific weight distributions P{w). We begin with the 
well-studied exponential disorder function f{x) = e^^ , where a; is a random number between 
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and 1 . From Eq. ^ follows that P{w) = l/{aw), where w E [l,e°-]. Using Eq. ^ 
we obtain Wc = e"'^", and from Eq. (jSJ, A becomes 

A = apc (12) 

For fixed S, we expect to obtain the same optimal path behavior. Indeed, this has been 



shown to be valid 



19|. 



Next we study A for the disorder function f{x)= x"", with x between and 1, where a 



is a parameter which can be positive or negative 



. For this case the disorder distribution 



is a power law P{w) = \a\ ^. Following Eqs. (jH]) and 0, we obtain 

A = \a\ (13) 

Note that a plays a similar role to a in Eq. ()12j) . but here S is independent of pc, which 
means that networks with different pc, such as ER networks with different average degree 
< k >= 1/pc, will yield the same optimal path behavior. 

We further generalize the power law distribution with the disorder function f{x) = x"' by 
introducing the parameter < A < 1 which is defined as the lower bound of the uniformly 
distributed random number x, i.e., 1 — A < a; < 1. Under this condition, the distribution 
becomes 

PM = ^^ ^G[(l-A)M]. (14) 
Again using Eqs. (jHI) and (jni),we obtain 

A= , (15) 

Apc + 1 - A ^ ^ 

Table H] shows the results of similar analyses for the lognormal, Gaussian, uniform and 
exponential distributions P{w). From Table HI we see that for exponential function, power 
law and lognormal distributions, A and therefore S can become large. However for uniform, 
Gaussian and exponential distributions, A is limited to a value of order 1, so S ^ 1 for 
large L. Note that for these distributions, A is independent of a or a. Thus only the former 
distributions (exponential function, power law, lognormal) are broad enough, and can lead 
to strong disorder behavior for the optimal path, while for the latter distributions (uniform, 
Gaussian and exponential), 5* is always small and only weak disorder can appear. 

To test the validity of our theory, we perform simulations of optimal paths in 2d square 
lattices and ER networks. Random weights from different disorder functions were assigned 



to the bonds. For an. L x L square lattice, we calculate the average length i of the optimal 
path from one lattice edge to the opposite. For an ER network of nodes, we calculate i 
between two randomly selected nodes. 

Figure 1 shows our simulation results for a 2d square lattice for power law and lognormal 
distributions. For small L, l/L^''^'^ approachs a constant, indicating the strong disorder 
regime ^ ~ L^'^^, while for large L it approaches to indicating i ^ L. The scaled 

results shown in Fig. lb show collapse for the same S for all distributions and therefore 
support Eq. (jH)) that S controls the optimal path, and that the optimal path i behaves in 
the same universal way for any broad distribution. Note that for large S~'^ the scaled curves 
approach a slope of 1 — dopt — —0.22 indicating the range of weak disorder. 

Scaled distributions p{i) for general disorder distributions P{w) are shown in Fig. 2, in 
which L, a, A and a are selected to keep 5* constant. Figure 2 suggests that for the same 
value of S, the distributions p{£) collapse when plotted versus the scaling parameter i/A'"^°pK 

Simulations for optimal paths on ER networks are shown in Fig. 3. Here we use the 
bombing algorithm P, to determine the path length ioo in the strong disorder limit, 
which is related to by £oo ~ iV"*" = N^^^ 0,13 • We see that for all disorder distributions 
studied, i scales in the same universal way with S^-^ = ioo/^- For 5^1, i/A is linear with 
log(^ooM) as expected (Fig. 3a). For large S = A/i^o (Fig. 3b), £ oc £oo ~ N'^^^, which is the 
strong disorder behavior 7|. Thus, we see that when increases, a crossover from strong 
to weak disorder occurs in the scaled optimal paths i/A vs. 5*^^. Again, the collapse of all 
curves for different disorder distributions for ER networks supports the general condition of 
Eq. dHD. 

Next we use Eq. (jH} to analyze the other types of disorder given in Table HI that do 
not have strong disorder behavior. For a uniform distribution, P{w) = 1/a and we obtain 
A = 1. The parameter a cancels, so = L~^/'^ for lattices, and S = N~^^^ for ER networks. 
Hence for any value of a, S" <^ 1, and strong disorder behavior cannot occur for a uniform 
distribution. 

Next we analyze the Gaussian distribution. We assume that all the weights Wj are positive 
and thus we consider only the positive regime of the distribution. We obtain A as given in 
Table I. The disorder is controlled solely by Pc which is related only to the type of network, 
and A cannot obtain large values. Thus, also for the Gaussian P{w), all optimal paths are 
in the weak disorder regime. 
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Similar considerations lead to the same conclusion for the exponential distribution. Sim- 
ulation results for the Gaussian distribution shown in Fig. 4 display only weak disorder (i.e. 
no weak-strong disorder crossover), thus supporting the above conclusions. 

In summary, we find a criterion for the disorder strength S on the optimal path in weighted 
networks for general distributions P{w). We find an analytical expression, Eq. (jS)), which 
fully characterizes the behavior of the optimal path. We simulated several distributions and 
the results support our analytical prediction. It is plausible that the criterion of Eq. (jH)) 
is valid also for other properties in weighted networks — such as conductivity and flow in 
random resistor networks — due to a recently-found close relation between the optimal path 
and flow 

Our results also suggest the conjecture that disorder distributions fall into two classes: 
those P{w) that can lead to strong disorder and those that cannot. Our studies suggest 
that the ratio cr//i determines the class, where o"^ is the variance of P{w) and fi the mean 
value. If a/fi can become large 1), strong disorder can occur, but if cr/yU ^ 1 only weak 
disorder occurs. This criterion is valid for all the cases simulated and we conjecture that it 
is valid for any weighted distribution. 

We thank ONR, European NEST project DYSONET, and Israel Science Foundation for 
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TABLE I: Parameters controlling the optimal paths on networks for various distributions of disor- 
der. 



Name 


Function 


Distribution 


A 


Livcrsc 




ir e [l.f") 

au; L / 




Power Law 




^1/a-l 

|a|A 


Aapc 
Apc+l-A 


Lognormal 


g\/2(Terf-i(2a;-l) 


g-(lniL02/2<T2 


\/2ttPc<j 




e-[erf-i(2pc-l)]^ 


Uniform 


ax 


1/a 


1 


Gaussian 


^f2a&xi-^{x) 


2e-ti-2/(2<T2) 






2erf-i(pc) 


Exponential 


C 




(pc-l)ln(l-pc) 
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FIG. 1: Simulation results for different disorder distributions in 2d lattice (a) before and (6) after 
scaling. A lower limit of L = 15 is set to avoid the finite size effects, (a) Curves for power law 
disorder chosen from x"", < x < 1, with a from 20 to 40, fixing a = 40 and A = 0.6 or 0.8, 
and lognormal distribution with a from 5 to 20. (6) Same results plotted in a scaled form. Two 
lines are drawn in (6), with slopes and —0.22, indicating the theoretical behavior of optimal path 
in the strong and weak disorder limits respectively. The inset in (b) shows the ratio between the 
highest weight and total weight versus S^"^ . A line of Wmax/wtotai = 0.5 indicates the point S* 
above which the highest weight dominates the total weight. 
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l/A °p' 

FIG. 2: Scaled distributions of p(^) in 2d lattice. Distributions used here are lognormal distribution 
with G = 20, power law distribution with a = 30 and power law distribution with a = 40, A = 0.6. 
The linear size L of lattice is selected by fixing S = 0.8. 
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FIG. 3: The function i/A for ER networks after scaling, where (a) is a hnear-log plot and {h) 
is a log-log plot. Distributions used are power law x°' with 10 < a < 30 where < x < 1, re" 
with a = 30 and the range of A < j; < 1 with A = 0.6 or 0.8, and lognormal distribution with 
10 < cr < 30. The straight line in (a) indicates weak disorder and the straight line in (6) indicates 
strong disorder. 
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FIG. 4: The optimal path for Gaussian distribution of weights, (a) Gaussian distribution for ER 
networks and (6) Gaussian distribution for 2d lattices. Note that these curves would collapse after 
scaling to the curves in figure lb in the weak disorder tail of large S~'^ . 
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